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Abstract 

Finite temperature phase boundary between superfluid phase and nor- 
mal state is analytically derived by studying the stability of normal state 
in rotating bosonic optical lattice. We also analytically prove that the 
oscillation behavior of critical hopping matrix directly follows the upper 
boundary of Hofstadter butterfly as the function of eff'ective magnetic 
held. 



PACS number(s): 03.75.Lm, 05.30.Jp, 73.43.Nq 



1 Introduction 

Bose-Hubbard model of interacting bosons on a lattice has been used to describe superfluid- 
Mott insulator (MI) phase transition in a variety of systems at zero temperature, e.g., 
Josephson arrays and granular superconductors [Ij. The recent suggestion to experimen- 
tally observe this transition in a system of cold bosonic atoms in an optical lattice |2j and 
its successful experimental demonstration [3] have aroused much theoretical [H IH E] and 
experimental [TJ |8] interest in this model, especially rotating optical lattice has also become 
brand-new topics in bosonic system. Most of work about rotating optical lattice focused 
on the superfluid phase and studied the pinning effect of the vortex lattice due to optical 
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periodic potential [9l [ini [HI [E] . The similar question has been investigated in the type-II 
superconductor [13l [U [ISl [H] . 

However the question of superfluid-MI phase transition at zero temperature still exists 
in rotating optical lattice and the phase diagram at zero temperature has been achieved by 
strong coupling expansion [T7] and Gutzwiller approach [THl [TO] ■ Strong coupling expan- 
sion obtained the phase boundary between the superfluid phase and MI by perturbatively 
computing the energy difference between MI and single-charge excitation states on top of 
the MI. This method is very accurate and applicable to the random dimension but is not 
suitable for system in deep superfluid phase where the perturbation is not valid any longer. 
Gutzwiller approach is at the self-consistent mean-field level and based on an ansatz that 
many-body ground state factorizes into product of single lattice site wave function. So under 
this approximation the system become diagonal with respect to the lattice site and we can 
use an effective single-site Hamiltonian. The disadvantage of Gutzwiller approach is that it 
fails to describe the correct short-range correlation between different lattice sites and so is 
an uncontrolled approximation. However, Gutzwiller approach can predict a qualitatively 
similar phase diagram with strong coupling expansion [18]. 

In this paper, we mainly extend the phase diagram of a rotating two-dimensional bosonic 
optical lattice to finite temperature utilizing the Gutzwiller mean-field theory. At finite 
temperature, MI is replaced by normal state which possesses finite compressibility. Here 
we do not include the crossover from MI to normal state at finite temperature since there 
is not an conventional definition for this crossover as far as what we know is concerned. So 
at finite temperature the phase transition happens between superfluid phase and normal 
state instead of MI. In section 2, by making an analogy between rotating optical lattice and 
electrons constrained by periodic potential and external magnetic field, we qualitatively 
derive the Bose-Hubbard Hamiltonian in rotating optical lattice. In section 3, we follow 
the method used in [20] to analytically locate the phase boundary of superfluid phase and 
normal state by discussing the stability of flxed point corresponding to the normal state 
and in section 4 a brief conclusion is given. 

2 Bose-Hubbard Hamiltonian in Rotating Optical Lattice 

We consider a two-dimensional bosonic system in XY plane restricted by a square optical 
lattice and an harmonic trapping potential which have a common rotating velocity along 
Z axis. In the laboratory frame, the potential is generally time-dependent. It is therefore 
convenient to transform to the frame rotating with the potential, since in that frame the 
potential is constant in time, and thus the standard methods for finding the equilibrium 
may be employed. In the frame of rotating potential, the second quantized Hamiltonian 
for a particle of mass m in an harmonic trap of natural frequency u; can be written into 
H = Ho + Hi with 



Ho 



I 



( 




(1) 




(2) 



2 



where Lz = —ifi{XdY — Ydx) is angular momentum along Z axis and g{> 0) is the strength 
of contact interaction potential between two particles. Vl{R) is periodic optical potential 
and ^'(-R) is field operator for boson particles. The term involving the chemical potential 
^ is added because it is very convenient to be in the grand canonical ensemble. Hq can be 
rearranged into 



-ihV - e/cA{R)f 
2m 



+ VL{R) + ^m{io^ 



n^)R'^ - fi 



^{R) (3) 



where A{R) = mc/eO. x i? is an effective vector potential with c and e representing light 
speed and charge quanta. This form suggests that the effects of rotation are partitioned 
into two different parts. The term l/2m(a;^ — Q?')R? implies that the centrifugal potential 
weakens the role of trapping potential (in order to stabilize the system, Vt <lo). The other 
part of rotation is included in the first term whose role is producing an effective magnetic 
field i? = V X A{R) = 2mc/eQ and provides a structure of Landau level for atoms. Hence 
at this point we can draw a conclusion that the motion of atoms in a rotating optical lattice 
under the assumption w = is completely the same as that of electrons constrained by 
periodic potential and external magnetic field. Therefore the method utilized to deal with 
electrons can be applicable to the our problem. For simplicity, below we assume lo = 0, so 
that the centrifugal force accurately compensates the harmonic trapping potential. 

For systems without rotation and trap potential (0 = a; = 0) we usually assume that 
the atoms are cooled within the lowest Bloch band and the field operator can be expanded 
into ^{R) = ^ihiw{R — Ri) in terms of the lowest Wannier function w{R — Ri) with bi 
being bosonic annihilation operator at the lattice site Ri. Hence the system is sufficiently 
described by a single-band Bose-Hubbard Hamiltonian [2] 



(4) 



<ij> 



where t = -JdRw*{R- Ri) 
to nearest neighbors and U = 



+ Vl{R) 'w{R — Rj) is hopping matrix restricted 
Ri)\'^ is on-site interaction strength. When 



2m 

g J dR \w{R — ±vi) 

external effective magnetic field appears, many-band effects must be considered which leads 
to that single-band Bose-Hubbard Hamiltonian is not valid any more. Fortunately from the 
study of lattice electron in external magnetic field [21] we know that there exists an effective 
single-band Hamiltonian which can be obtained from @ by only making a substitution for 
hopping matrix 



H 



-t exp 

<ij> 



-- / A{R) ■ dR 
he Jr. 



U 



(5) 



The added phase factor is called Peierls phase factor. It is well known that superfluid-MI 
phase transition comes from the competition between hopping matrix and on-site interaction 
strength [3]- But when the magnetic field is introduced the on-site interaction strength is 
unchanged while the hopping matrix is modified, so we naturally expect a modified phase 
boundary between superfluid and MI. 
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The vector potential A[R) in the symmetric gauge concerns the X and Y components of 
coordinate at the same time and makes below calculation be more complex. We hope only 
X or Y component is concerned, which is realized by making a canonical transformation to 
field operator 



bi exp 



he 



Ri 



mc/efl X R ■ dR 



(6) 



with R = {X, —Y) and Ri being a arbitrary reference point. Under this transformation 



H = —t ^ exp 

<i3> 



i27r(/) / xdy 



bib, 



+ 



^T^blblbA 



(7) 



where all the coordinates are scaled by the lattice constant a hence are dimensionless. 
(j) = Ba^ / {he/ e) represents the number of magnetic flux quanta penetrating the unit cell. 
In fact it is easy to find that the above Hamiltonian corresponds to that in Landau gauge 
A{R) = BXY with Y denoting the unit vector along the Y axis. In the next section, we 
will regard the Hamiltonian ^ as our starting point and study its phase diagram at finite 
temperature at mean-field level. 



3 Phase Diagram in Gutzwiller Mean-Field Approach 

The Gutzwiller approach is a self-consistent mean-field method and equivalent to the de- 
coupling approximation |22l |23] to the hopping term 



b\hj 



< b\ > bj + bl <bj > - < bl X bj > 



Oibj + b-ttj 



(8) 



where a,- 



a. 



is superfluid order parameter which distinguishes the superfluid phase from 



normal state. If magnetic field vanishes, the whole system is uniform and order parameter 
is also site-independent. But we can not suppose this when the magnetic field appears. 
After this decoupling, the system is describable in terms of a single-site Hamiltonian 



-t 



U 



blimbnmbnmbnm ~ fJ'blimbnm 



■i2TTn<j> 



n(m+l) + e*^''"''^«n(m-l) ) + H.C. 



+ 



(9) 



where we label site of the lattice i by two ordered integers i = (n, m), the first integer along 
the X axis and the second one along the Y axis. In the Landau gauge, hopping along the 
Y axis achieves the Peierls phase factor and that along X axis is invariant. The above 
Hamiltonian has two striking characteristics [18j . On the one hand, it is independent of Y 
component, so the translational symmetry along Y axis is conservative and we can suppose 
order parameter anm = ct„ in correspondence with the case without magnetic field. On 
the other hand although it depends on X component, for rational 4> = p/q {p,q have no 
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common factor), q-site translational symmetry along X axis is recovered = ctn+q- So 
the Hamiltonian is further reduced into 

Hn = -t 



bl (a„+i + an-i + 2a„ cos 27rn(/)) + H.C. + —hlhlhnhn - fJ-bibn (10) 

with n being integer from 1 to q. In addition, the Hamiltonian is periodic as the function 
of magnetic field Hn{(l)) = Hn{4> + K) with K being a random integer so that we only need 
consider (p £ [0, 1). Note also that in Q and (|10|) we have neglected a constant term which 
does not influence our result. 

The self-consistency of Gutzwiller method must be carried out by the condition 



a-n = -^Tr (bn 



fill 



with /3 = 1/{KbT) and partition function Z„ = Tr exp {—l3Hn)- Introducing the same 
notation in ^20j 7„ = tan, self-consistent condition can be rewritten into 



1n = 



d 



+ 



d 



+ 



1 



d 



6/3 \(?7n+i djn-i 2 COS 27rn(/) (?7; 



In Zr, 



(12) 



Under Gutzwiller approximation the eigenstates of Hn can be expanded in terms of Fock 
state, then we diagonalize this Hamiltonian matrix under Fock basis truncated until a given 
number of particle N to obtain the eigenvalues Ei {I = 0,1,- ■ ■, N) which are function of 
7n-i,7n and 7n+i- So we proceed using the eigenvalues 



In 
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N 
1=0 



dEi dEi 

+ - — - + 



1 



djn-i d^in+i 2cos27rn097, 



dEA ^ 



(13) 



The derivative of eigenvalue is computable from the characteristic polynomial P{\, 7„_i, 7„, 7n+i) 
of Hamiltonian matrix |24j 



9P(£';,7n-l,7n,7n+l) 

OEi 

<9-P(£^Z,7n-l,7n,7n+l) 



dP{Ei,-fn-l,7n, 7n+l) ,dP{Ei , Jn-1, In, In+l) 



I- 



dEr 



'n-\ 



N 

Y,P^^HEi,-fn-l,ln,ln+l) 
k=0 

2(7n-i + 7n+i + 27„ COS 27rn<; 
N 



-l,ln,ln+l) 



lk=l 



dP{Ei,Jn-l,ln,7n+l) 

d-fn+1 
dP{Ei,-^n^l,ln,ln+l) 
djn 



dP{Ei,Jn-l,7n,7n+l] 



2 cos 2Tm<f) 



dP{Ei,-fn-l,ln,ln+l) 



n-l 



(14) 



where the notation P^^"^ {Ei,'yn-i,ln,ln+i) denotes the characteristic polynomial of the 
matrix obtained by discarding from Hamiltonian matrix the rows and columns labeled 
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by the set of indices {k} and the polynomial (5(7n-i5 7ni 7n+i) satisfies (5(0,0,0) = 0. 
Substituting all above relations into (jl3p 



7n = H7n-1 + 7n+l + ^In COS 27rn0) 



E^=l A:P(^-''-lH^«,7n--l,7n,7n+l) + Q(7n-1. 



VF(^;,7„_l,7„,7„+l) = ^jV r,(k)CP ^ 

Efc=0 7n-l, 7n, 7n+l) 

Until now, we have obtained the equation set of self-consistently deciding all order param- 
eters. If order parameters have nonzero solution the system is in superfluid phase. If order 
parameters have zero solution, the system is in normal state. So we can determine the phase 
boundary between superfluid phase and normal state by studying the stability of the flxed 
point corresponding to the normal state (7^ = for all n) [20]. According to a standard 
theory, the stability of such a fixed point can be discussed based on the spectrum of the 
matrix linearizing the map defined by equation set (llSp in the vicinity of normal state. 
Linearizing (jlSp around the fixed point of normal state, we obtain 

7n = t6([/,;U,/3)(7„_i +27„cos27rn(/) + 7„+i) (16) 
a,,jr _ Y.loe-^^^''^^^-^)-^^)W{U/2l{l - 1) - 0,0,0) 



Writing above equations in the form of matrix if we denote T = (71,72, • • -75 



M{4>) 



te{U,fi,l3)M{(f))T 

( 2 cos 27r(/> 1 • • • 1 

1 2 cos 47r(/) 1 • • • 

1 • • • 1 2 cos 27rg(/> j 



(18) 



Therefore, the fixed point of the normal state is stable only if the maximal eigenvalue 
£maxi4') of M((l)) satisfies tQ{U, fi, f3)emax{4') < 1) which signifies that the phase transition 
between superfluid phase and normal state is 

teiU, fl, (3)€maA<P) = (19) 

This is our main result. At T = OK only the eigenstate having the lowest energy contributes 
to the partition function. Easily proven when // G [Z — 1, /], the phase boundary at T = OK 
is specifled by 

emax{(t>) U + fl 

which is the same as the result in \W\ I25| . 

Below we connect emax{4>) with the famous Hofstadter butterfly. We find that eigen- 
values of the matrix M{(f)) actually correspond to a part of energy spectrum of Hofstadter 
butterfiy [26j. According to the proof in [19j that the maximal eigenvalue of M{(j)) is equal 
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to the maximal eigenvalue of Hofstadter butterfly, hence from (I19p we find for fixed U, /i, /? 
the critical hopping strength is inversely proportional to the maximal eigenvalue of Hofs- 
tadter butterfly. According to the fact that the maximal eigenvalue of Hofstadter butterfly 
shows an oscillatory behavior as the function of magnetic field [26], the critical hopping 
strength also exhibits the oscillation following the maximal eigenvalue, i.e. upper boundary 
of Hofstadter butterfly. 




0.0 



Figure 1: The three-dimensional phase diagram for different temperature KbT/U = 
0.0, 0.03, 0.08. Below the critical surface is normal state. 

In Fig.l, we plot three-dimensional phase diagram for different temperature. When 
T = OK the phase boundary separates the superfluid phase and MI, which is not true for 
finite temperature. In order to concentrate on the effects of temperature and magnetic 
field, we also plot the two dimension phase diagram for different temperature and magnetic 
field in Fig. 2. Note that the phase boundary of = is obtained by letting emax(O) = 4 
[25\ I26j. Seeing from the Fig. 2, we could draw below conclusions. For fixed magnetic field 
the higher the temperature, the smaller the area of the superfluid phase, which is consistent 
with the fact that the high temperature destroys the superfluidity. At the same time for 
fixed temperature, magnetic field has a much more complex effect on the superfluidity in 
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view of oscillatory behavior of critical hopping matrix. Generally speaking in contrast to the 
situation without the magnetic field, magnetic field always nonmonotonically increases the 
area of normal state, which can be illustrated from the bandwidth of Hofstadter butterfly. 
On the one hand the narrower the bandwidth, the smaller the effective hopping matrix. 
But the bandwidth of Hofstadter butterfly is often less than that without magnetic field 
[26| . Hence the effective hopping matrix is always less than t. On the other hand the 
on-site interaction strength is invariant. Considering the above two factors, we naturally 
understand the effect of magnetic field. 

4 Conclusion 

In conclusion we have qualitatively derived the Hamiltonian of rotating optical lattice, ana- 
lytically extended the phase diagram of rotating Bose-Hubbard model to finite temperature 
and analyzed the relation between the oscillation behavior of critical hopping matrix and 
Hofstadter butterfly. In addition, we have also illustrated how the rotation influences the 
phase diagram. 
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Figure 3: The three-dimensional phase diagram for different temperature KbT/U = 
0.0, 0.03, 0.08. Below the critical surface is normal state. 



Figure 4: The two-dimensional phase diagram for different temperature and magnetic field. 
Below the critical line is normal state. 
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